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The  performances  of  thermal-concentrated  solar  thermoelectric  generators  (TEGs)  at  three  different 
geometric  types  are  investigated  numerically  to  aid  in  designing  practical  devices.  The  temperature- 
dependent  properties  of  the  commercial  thermoelectric  material  are  taken  into  account,  and  an  equiv¬ 
alent  model  based  on  the  three-dimensional  finite  element  scheme  is  developed  to  simplify  and 
accelerate  simulations.  The  constriction  thermal  resistance  and  thermal  spreading  resistance  are 
considered  in  the  equivalent  model.  Increasing  substrate  area  increases  the  thermal  concentration  ratio; 
this  improves  the  performance  of  the  solar  TEG.  In  the  three  geometric  types,  the  smallest  element  with 
the  substrate  area  of  90  x  90  mm2  provides  the  maximum  system  efficiency  of  4.15%.  For  a  TEG  at  a  given 
element  length,  decreasing  the  cross-sectional  area  of  the  thermoelectric  element  is  a  feasible  route  to 
improve  the  performance.  Under  the  situation  of  forced  convection,  varying  convection  heat  transfer 
coefficient  has  an  ignorable  effect  on  the  performance.  For  equal  convective  heat  transfer  coefficients, 
water  cooling  is  better  than  air  cooling  for  the  net  output  power  of  the  TEG  because  of  its  increased 
specific  heat.  Therefore,  water  cooling  is  recommended  for  the  cooling  of  the  solar  TEG. 

©  2013  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Thermoelectric  generators  (TEGs)  are  energy  conversion  de¬ 
vices,  which  unlike  convectional  heat  engines,  are  entirely  solid- 
state,  environmentally  friendly,  extremely  reliable,  simple, 
compact,  and  safe  [  1,2],  When  a  temperature  difference  across  a  TEG 
exists,  it  will  generate  electrical  power  through  the  Seebeck  effect 
[3,4],  Although  the  principle  of  thermoelectricity  was  discovered  in 
1834,  there  were  few  practical  applications  until  the  middle  1950s 
[5],  Since  the  first  and  the  second  energy  crises  occurred  in  the 
1970s,  considerable  interest  has  arisen  in  the  development  of 
thermoelectric  power  generation  [6],  The  applications  of  TEGs  are 
limited  due  to  their  relatively  low  conversion  efficiency  [7,8],  Waste 
heat  is  extensively  encountered  in  many  industrial  and 
manufacturing  processes  [9,10]  and  it  is  low-cost  and  even  no-cost 
if  waste  heat  can  be  recovered.  Therefore,  the  low  efficiency  prob¬ 
lem  of  TEGs  is  not  a  critical  issue  when  they  are  applied  in  such 
settings  [7]  if  waste  heat  is  used  as  the  heat  source  of  the  TEG. 
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Moreover,  the  energy  conversion  process  of  TEGs  is  totally  green 
[11  ]  so  that  they  are  able  to  abate  greenhouse  gas  emissions. 

A  number  of  practical  applications  of  thermoelectric  devices 
through  waste  heat  recovery  can  be  found  in  literature.  Hsiao  et  al.  [  1 2  ] 
proposed  a  theoretical  model  to  predict  the  performance  of  a  TEG 
module  by  recovering  waste  heat  from  an  automobile  engine.  Their 
results  showed  that  the  performance  of  the  TEG  installed  on  an 
exhaust  pipe  was  better  than  that  installed  on  a  radiator.  Hsu  et  al.  [  1 3  ] 
constructed  a  thermoelectric  system,  which  comprised  24  TEG 
modules  to  recover  waste  heat  from  the  exhaust  pipe  of  an  automo¬ 
bile,  to  find  its  optimum  operating  conditions.  Chen  et  al.  [14]  inves¬ 
tigated  the  characteristics  of  thermoelectric  modules  used  for  power 
generation  by  recovering  low-temperature  waste  heat  at  various 
operating  conditions.  They  discovered  that  the  effects  of  the  flow 
pattern  of  the  heat  sink  and  the  water  flow  rate  on  the  performance  of 
the  modules  were  not  significant,  whereas  the  heat  source  or  the 
heating  temperature  played  an  important  role  on  the  performance. 

Another  low-cost  heat  source  is  the  solar  energy  which  has  been 
widely  applied  in  industries,  such  as  solar  thermal  systems  and 
photovoltaic  systems  [15],  In  recent  years,  solar-driven  TEGs  have 
received  a  great  deal  of  attention  and  a  lot  of  studies  have  been 
carried  out.  One  of  the  common  designs  of  solar  thermoelectric 
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Nomenclature 

Ql 

heat  transfer  rate  at  the  solar  TEG’s  cold  side  (W) 

QL 

thermal  load  vector  (W) 

Ac 

cross-sectional  area  of  the  collector  (mm2) 

T 

temperature  (°C) 

Cp 

specific  heat  at  constant  pressure  (kj  kg-1  K-1) 

Te 

vector  of  nodal  temperatures  (°C) 

Copt 

optical  concentration  ratio  of  solar  thermoelectric 

Tf 

surface  temperature  of  fins  (°C) 

generator 

tf 

fin  thickness  (mm) 

Oh 

thermal  concentration  ratio  of  solar  thermoelectric 

Too 

environment  temperature  (°C) 

generator 

U 

mean  velocity  of  fluid  (m  s_1) 

D 

depth  of  thermoelectric  element  (mm) 

W 

width  of  thermoelectric  element  (mm) 

Dg 

fin-to-fin  spacing  (mm) 

WHS 

width  of  heat  sink  (mm) 

E 

electric  field  intensity  vector  (V  m  1 ) 

F 

area  factor 

Greek  letters 

C 

curve  fitting  functions  of  relative  error 

a 

Seebeck  coefficient  (V  K_1) 

Hf 

fin  height  (mm) 

T 

electrical  conductivity  (S  m  ') 

hL 

convection  heat  transfer  coefficient  at  the  TEG’s  cold 

£ 

emissivity  of  the  solar  absorber 

side  (W  m“2  K1) 

Pa 

absorptivity  of  the  collector  coating 

IL 

electric  current  load  vector  (A) 

Popt 

optical  efficiency  of  the  Fresnel  lens 

1 

electric  current  density  vector  (A  m-2) 

Vt 

total  efficiency  of  the  system 

ktt 

thermal  stiffness  matrix 

Vteg 

conversion  efficiency  of  the  TEG 

K'i’T 

Seebeck  stiffness  matrix 

A 

fluid  viscosity  (N  s  m-2) 

K<p<p 

electric  stiffness  matrix 

P 

fluid  density  (kg  m-3) 

k 

thermal  conductivity  (W  m  1  K-1) 

Pe 

electrical  resistivity  (Q  m) 

L 

length  of  thermoelectric  element  (mm) 

a2 

coefficient  of  determination 

Lhs 

length  of  heat  sink  (mm) 

OSB 

Stefan-Boltzmann  constant 

m 

mass  flow  rate  of  fluid  (kg  s  ') 

(  =  5.67  x  10-8Wm-2r4) 

m 

constant  for  the  curve  fitting  function 

0 

electric  scalar  potential  (V) 

N 

element  shape  function 

<Pe 

vector  of  nodal  electric  potentials  (V) 

rtf 

number  of  fins 

Ap 

pressure  drop  across  the  heat  sink  (N  m  2) 

n 

constant  for  the  curve  fitting  function 

^equivalent  output  power  of  equivalent  model  (W) 

Subscripts 

Pout 

output  power  (W) 

a 

air 

Epump 

pumping  power  (mW) 

f 

fin 

Ereal 

output  power  of  real  model  (W) 

H 

hot  side 

Qs 

solar  irradiance  (W  m-2) 

HS 

heat  sink 

Q 

heat  transfer  rate  (W) 

L 

cold  side 

Qh 

heat  transfer  rate  at  the  solar  TEG’s  hot  side  (W) 

rd 

radiation 

Qin 

input  energy  of  solar  thermoelectric  generator  (W) 

w 

water 

systems  is  the  optical-concentrated  solar  TEG  which  uses  an  optical 
concentration  system  to  focus  sunlight  on  the  hot  side  of  the  TEG.  Li 
et  al.  [16]  designed  a  prototype  optical-concentrated  solar  TEG  and 
used  a  numerical  method  to  evaluate  its  performance.  Their  results 
showed  that  the  efficiency  of  the  TEG  could  reach  up  to  9.8, 13.5, 
and  14.1%  when  Bi2Te3,  skutterudite,  and  LAST  alloys  were 
employed  as  the  materials  of  the  TEG,  respectively.  Fan  et  al.  [17] 
fabricated  and  experimentally  tested  an  optical-concentrated  so¬ 
lar  TEG  system  along  with  a  parabolic  dish  concentrator.  The  sys¬ 
tem  was  able  to  produce  electric  power  up  to  5.9  W  under  the 
temperature  difference  of  35  °C  with  the  hot-side  temperature  of 
68  °C.  Yang  et  al.  [18]  developed  a  numerical  model  to  analyze  the 
performance  of  a  TEG  module  with  optical  concentration.  Their 
results  suggested  that  the  efficiency  of  the  TEG  decreased  to  3.85% 
and  had  a  relative  drop  of  61.3%  when  the  contact  resistance  and  all 
heat  losses  were  considered.  Xiao  et  al.  [19]  analyzed  the  perfor¬ 
mance  of  a  multi-stage  solar  TEG  by  means  of  a  numerical  method. 
They  pointed  out  that  that  the  efficiency  of  the  solar  TEG  system 
could  reach  10.52%  based  on  a  three-stage  thermoelectric  module. 

As  far  as  the  thermal-concentrated  solar  TEG  is  concerned,  a 
flat-panel  solar  absorber  on  the  hot  side  of  the  TEG  is  used  to 
concentrate  the  solar  heat  through  heat  conduction.  Kraemer  et  al. 
[20]  designed  a  flat-panel  solar  TEG  system  with  high  thermal 
concentration;  the  system  with  the  conditions  of  a  thermal  con¬ 
centration  ratio  of  299  and  AM1.5G  (1  kW  m-2)  could  attain  a  peak 


efficiency  of  4.6%.  Chen  [21]  developed  a  theoretical  model  to 
evaluate  the  efficiency  of  a  TEG  pair  in  association  with  thermal 
concentration  and  optical  concentration.  His  results  suggested  that 
the  efficiency  of  the  TEG  pair  under  thermal  concentration  but  with 
little  or  no  optical  concentration  could  be  larger  than  5%  when  the 
hot-side  temperature  was  between  150  and  250  °C. 

According  to  the  above  literature  review,  the  two  routes  of  op¬ 
tical  concentration  and  thermal  concentration  can  be  utilized  to 
improve  the  performance  of  solar  TEGs.  The  optical  concentrator  is 
generally  costly  [22],  On  the  other  hand,  very  little  research  has 
been  performed  on  the  thermal-concentrated  solar  TEG,  especially 
in  numerical  simulation.  For  this  reason,  the  objective  of  this  study 
is  to  numerically  investigate  the  performance  of  a  solar  TEG  with 
thermal  concentration.  A  three-dimensional  finite  element  scheme 
will  be  employed  to  model  the  TEG  system.  In  some  studies 
[12,13,21  ],  the  properties  of  thermoelectric  materials  were  assumed 
to  be  constant.  This  may  induce  a  numerical  deviation  at  high 
temperatures.  Consequently,  the  temperature-dependent  proper¬ 
ties  of  the  thermoelectric  materials  are  considered  in  this  work.  In 
particular,  an  equivalent  model  is  developed  to  simplify  and 
accelerate  the  numerical  approach.  The  effects  of  substrate  area  and 
thermoelectric  element  geometry  on  the  performance  of  the  TEG 
will  be  explored.  Eventually,  the  power  consumption  for  the  cool¬ 
ing  of  the  TEG  is  predicted  via  a  theoretical  method  to  evaluate  the 
net  output  power  of  the  TEG  system. 
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2.  Methodology 

2.1.  Physical  model  and  assumptions 

The  schematic  of  the  investigated  thermal-concentrated  solar 
TEG  system  is  shown  in  Fig.  1  a.  In  the  system,  a  solar  absorber  with 
the  heat  uptake  of  Qh  is  located  at  the  hot  side  of  the  TEG.  Mean¬ 
while,  the  heat  Ql  is  liberated  to  the  heat  sink  at  the  cold  side  of  the 
TEG.  The  generated  electrical  current  is  used  to  power  the  external 
loading.  Fig.  lb  demonstrates  the  practical  unpackaged  thermo¬ 
electric  module  with  the  size  of  100  x  100  mm,  which  contains  16 
thermoelectric  (TE)  pairs.  The  commercial  bismuth  telluride  ele¬ 
ments  were  purchased  from  a  company  in  Taiwan  (Tande  Energy 
and  Temperature  Associates  PTY.  CO.).  Prior  to  performing  simu¬ 
lations,  the  following  assumptions  are  adopted  to  simplify  the 
problem. 

(1)  The  solar  TEG  is  in  steady-state. 

(2)  The  configurations  of  the  p-type  and  n-type  elements  are 
identical. 

(3)  The  thermoelectric  elements  are  connected  electrically  in 
series  and  thermally  in  parallel. 

(4)  The  material  properties  of  the  thermoelectric  elements  are 
temperature-dependent,  but  the  material  properties  of  the 
solar  absorber,  substrate,  and  electrode  are  constants. 

(5)  The  whole  system  is  in  a  vacuum  environment  [20,21  ],  so  the 
heat  loss  due  to  heat  convection  is  neglected. 

(6)  The  radiation  heat  transfer  only  exists  between  the  hot-side 
surface  and  the  environment. 

(7)  The  solar  absorber  and  the  heat  sink  are  not  included  in  the 
computational  domain,  but  their  behavior  is  described  by 
boundary  conditions. 


Based  on  the  preceding  assumptions,  only  one  thermoelectric 
pair,  including  a  p-type  element  and  an  n-type  element,  is 
considered  in  the  computational  domain,  as  shown  in  Fig.  2a.  The 
substrate  area  for  thermal  concentration  is  usually  much  larger 
than  the  cross-sectional  area  of  the  elements.  For  example,  when 
16  TE  pairs  are  included  in  a  solar  TEG  module  with  the  substrate 
area  of  100  x  100  mm2,  the  average  substrate  area  for  each  TE  pair 
is  25  x  25  mm2,  whereas  the  cross-sectional  area  of  the  TE 
element  is  1  x  1  mm2,  as  shown  in  Fig.  2b.  The  thermal  concen¬ 
tration  is  characterized  by  the  thermal  concentration  ratio  (Cth) 
which  is  defined  as  the  area  ratio  of  the  solar  absorber  to  the 
cross-sectional  area  of  the  thermoelectric  elements  [20,21], 
Therefore,  the  value  of  Cth  in  Fig.  2b  is  312.5.  In  the  real  model,  the 
substrate  area  is  given  in  accordance  with  a  realistic  geometry.  In 
the  equivalent  model,  however,  the  substrate  area  is  assumed  to 
be  identical  to  the  given  element  geometry,  and  an  area  factor  is 
embedded  in  the  boundary  conditions  to  approach  the  input  en¬ 
ergy,  dissipated  energy,  and  the  radiation  heat  which  are  equal  to 
those  of  the  real  model.  The  advantage  of  the  equivalent  model  is 
that  the  number  of  computational  cells  can  be  minimized  so  as  to 
reduce  the  computational  time  significantly  when  compared  to 
the  real  model. 

2.2.  Governing  equations 

The  behavior  of  a  thermoelectric  material  is  described  by  the 
following  governing  equations  which  couple  the  thermal  and 
electrical  fields  [23,24]: 

V-(aTj'j -V-(kVT)  =  J-E  (1) 


(a) 


uuinui 

m  m  - 

rTTTTTTTTTTTTTTTTTl 


Fig.  1.  (a)  Schematic  of  thermal-concentrated  solar  TEG  system  and  (b)  photograph  of  the  unpackaged  thermoelectric  module  (P  and  N  represents  the  p-type  and  n-type  ther¬ 
moelectric  elements,  respectively). 
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(a) 


Real  model 


■  Electrode 


Equivalent  model 


(b) 


25  mm 

- 1 


25  mm 


Fig.  2.  Schematics  of  (a)  the  real  model  and  the  equivalent  model  and  (b)  the  top  view  of  a  TE  pair. 


(2) 

where  T  and  E  represent  the  temperature  and  the  electric  field, 
respectively.  The  electric jield  is  derived  from  the  electric  scalar 
potential  </>  ( E  =  -V0);  J  represents  the  electric  current  density 
vector;  a,  k,  and  pe  are  the  Seebeck  coefficient,  thermal  conduc¬ 
tivity,  and  electrical  resistivity  of  the  thermoelectric  element, 
respectively. 

2.3.  Boundary  conditions 

The  boundary  conditions  are  made  up  of  two  difference  sur¬ 
faces,  including  a  hot-side  surface  and  a  cold-side  surface. 

2.3.1  Hot-side  surface 

A  heat  flow  rate  Qjn  is  given  at  the  hot-side  and  it  is  obtained  by 

[16]: 

Qin  =  Qs  x  C0pt  x  Ac  x  7jopt  xrjaxF  (3) 

where  qs,  Copt,  Ac,  Jjopt,  Va.  and  Fare  the  solar  irradiance,  the  optical 
concentration  ratio,  the  substrate  area  per  TE  pair,  the  optical  ef¬ 
ficiency  of  the  Fresnel  lens,  the  absorptance  of  the  solar  absorber, 
and  the  area  factor,  respectively.  The  optical  concentration  ratio  is 
defined  as  the  area  of  the  lens  divided  by  that  of  the  solar  absorber 
[16],  The  cross-sectional  area  of  the  absorber  is  assumed  to  be 
equivalent  to  the  cross-sectional  area  of  the  thermoelectric  module. 
The  area  factor  only  exists  in  the  case  of  equivalent  model  and  it  is 
defined  as  the  substrate  area  per  TE  pair  of  the  real  model  divided 


by  that  of  the  equivalent  model.  Additionally,  the  heat  loss  from  the 
hot-side  surface  into  the  environment  is  described  by  heat  radia¬ 
tion  and  it  is  expressed  by 

Qrd  =  coSBAcF(T*  -  Ti)  (4) 

where  r,  <rSB,  TH,  and  T«,  are  the  emissivity  of  the  solar  absorber,  the 
Stefan— Boltzmann  constant  ( =  5.67  x  10-8  W  m-2  K-4),  the  hot- 
side  temperature,  and  the  ambient  temperature,  respectively.  In 
this  study,  the  material  of  the  solar  absorber  is  assumed  to  be 
aluminum  nitrogen  cermet  with  the  emissivity  (r)  of  0.08  and  the 
solar  absorptance  ( qa )  of  0.93. 

2.3.2.  Cold-side  surface 

Two  different  boundary  conditions  (Case  1  and  2)  at  the  cold- 
side  surface  are  considered.  In  Case  1,  a  constant  temperature  Tf 
is  given.  In  Case  2,  the  convective  heat  transfer  is  regarded  and  it  is 
expressed  as: 

-kVT  =  htF(T -Tm)  (5) 

The  convection  heat  transfer  coefficient  hL  is  assumed  to  be  uni¬ 
form  at  the  cold-side  surface.  It  is  noted  that  the  area  factor  F  is  only 
given  in  the  equivalent  model. 

2.4.  Numerical  method  and  system  performance 

The  governing  equations  are  discretized  by  a  finite  element 
scheme  based  on  the  Galerkin  method  [25],  Specifically,  the  phys¬ 
ical  scales  of  temperature  T  and  electric  scalar  potential  <p  over  a 
computational  element  are  approximated  by  [26]: 
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T  CO  T  CO 

Fig.  3.  Profiles  of  (a)  Seebeck  coefficient,  (b)  electrical  conductivity,  (c)  thermal  conductivity,  and  (d)  ZT  value  for  the  solar  TEG  at  Geometry  A. 


=  [N\{Te} 

(6) 

=  [N]{?e} 

(7) 

where  Te  ipe,  and  N  are  the  vector  of  nodal  temperature,  the  vector 
of  nodal  electrical  potentials,  and  the  element  shape  function, 
respectively.  By  integrating  Eqs.  (1)  and  (2)  based  on  the  Galerkin 
method,  the  differential  equations  can  be  written  as  algebraic 
equations  [26]: 

[S  £]{£}-{?} 

where  QL  and  IL  are  the  thermal  load  vector  and  the  electric  current 
load  vector,  respectively.  K17,  K<M',  and  KvT  are  the  thermal  stiffness 
matrix,  the  electric  stiffness  matrix,  and  the  Seebeck  stiffness  ma¬ 
trix,  respectively,  and  they  are  defined  as  [26]: 


K 17  =  J  VJV-[k]-VJVdV  (9) 


=  J  VN-^J-VNdV  (10) 


K<pt  =  J  vN-  j-JJ  •  |aj  -VN  dV  (11) 

The  three  matrices  and  two  load  vectors  are  solved  by  means  of  the 
commercial  software  ANSYS  vl2.0.1.  The  calculations  are  iterated 
until  the  convergence  criteria  are  satisfied.  Convergence  is  defined 
when  the  relative  differences  of  temperature  and  electric  potential 
between  two  iterations  are  smaller  than  10  6. 

After  the  simulation  is  completed,  two  indexes  are  employed  to 
evaluate  the  performance  of  the  system.  One  is  the  output  power 
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Pout  and  the  other  is  the  total  efficiency  of  the  system  rjT.  The  latter 
is  expressed  as  [18]: 

VT  =  tjTEG'^opft/a  (12) 

where  designates  the  conversion  efficiency  of  the  TEG  which  is 

obtained  from  ?jTEC  =  P0ut/Qm  [18].  The  material  properties  of  the 
thermoelectric  elements,  such  as  the  thermal  conductivity,  Seebeck 
coefficient,  and  electrical  resistivity,  are  measured  practically  and 
their  profiles  are  shown  in  Fig.  3.  The  values  of  the  thermal  con¬ 
ductivity  and  the  electrical  conductivity  of  substrate  and  electrode 
are  listed  in  Table  1.  The  contact  resistivity  is  set  to  be 
5  x  l<r10  O  m~2  [16], 


rhermoelectric  moduli 


* 


Thermoelectric  module 


2.5.  Analytical  model  of  power  consumption  of  the  cooling  system 


Fig.  4.  The  front  view  and  side  view  of  plate-fin  heat  sink. 


The  net  output  power  of  the  solar  TEG  is  evaluated  by  a  theo¬ 
retical  method.  The  configuration  of  the  plate-fin  heat  sink  is 
sketched  in  Fig.  4.  The  flow  direction  of  the  fluid  paralleling  the  heat 
sink  base  and  the  laminar  and  fully  developed  flow  is  assumed. 
Meanwhile,  the  fin  surfaces  are  postulated  to  be  isothermal  due  to 
the  simplicity  in  theoretical  analysis.  Therefore,  the  pressure  drop 
A p  across  the  heat  sink  is  obtained  by  [29]: 


A  P-2&H 


(13) 


where  p,  IHs,  Q,  and  Dg  are  the  fluid  viscosity,  the  length  of  the  heat 
sink,  the  mean  velocity  of  the  fluid,  and  the  fin-to-fin  spacing, 
respectively.  The  mean  velocity  is  derived  from  the  mass  flow  rate 
of  fluid  [29]  and  it  is  expressed  by: 


U  * 


m 

pDgHfn{ 


(14) 


at  the  conditions  of  qs  —  900  W  m-2, 7/opt  =  85%,  and  T„  =  25  °C 
where  only  one  thermoelectric  pair  was  simulated.  The  cross- 
sectional  area  and  the  length  of  the  element  were  1.3  x  1.3  mm2 
and  3.9  mm,  respectively.  Fig.  5a  suggests  that  the  current  pre¬ 
dictions  are  in  good  agreement  with  the  results  of  Li  et  al.,  thereby 
validating  the  numerical  method.  The  numerical  method  was  also 
compared  to  the  experimental  data  of  Kraemer  et  al.  [20]  where  the 
thermal-concentrated  solar  TEG  was  operated  at  the  conditions  of 
Cg  =  l,Cth  =  196,  J]opt  =  100%,  qs  =  1500  Wm~2,  and  T„  =  20°C, 
and  the  dimensions  of  the  thermoelectric  element  were 
1.35  x  1.35  x  1.65  mm3.  As  shown  in  Fig.  5b,  the  maximum  devi¬ 
ation  between  the  numerical  results  and  the  experimental  data  is 
10.5%.  The  deviation  may  be  attributed  to  the  adopted  material 
properties  from  literature  which  were  not  exactly  consistent  the 
real  materials  and  the  heat  loss  in  the  experiment  which  was  not 
considered  in  the  numerical  simulation.  The  above  comparisons 
reveal  that  the  developed  numerical  model  is  a  feasible  tool  to 
predict  the  performance  of  the  thermal-concentrated  solar  TEG. 


where  rh,  p,  Hf,  and  nf  are  the  mass  flow  rate,  fluid  density,  fin 
height,  and  the  number  of  fins,  respectively.  The  mass  flow  rate  is 
given  by: 


:P(Tf-T 


(15) 


where  Ql,  Cp,  Tf,  and  are  the  heat  dissipated  by  the  thermo¬ 
electric  module  at  cold  side,  the  specific  heat  at  constant  pressure, 
the  surface  temperature  of  fins,  and  the  environment  temperature, 
respectively.  Ql  in  each  simulation  is  solved  by  ANSYS.  Therefore, 
the  power  consumption  by  the  fan  or  pump  is  calculated  as  [29]: 


12  pLHSrh2 
P2DgHfn{ 


(16) 


3.  Results  and  discussion 

This  study  focuses  on  the  simulation  and  design  of  a  thermal¬ 
concentrated  solar  TEG.  Three  different  sizes  of  thermoelectric 
element  (W,  D,  L)  are  taken  into  account  for  comparison:  they  are 
Geometry  A  (1.0  mm,  1.0  mm,  1.6  mm),  Geometry  B  (1.4  mm, 
1.4  mm,  1.6  mm),  and  Geometry  C  (2.8  mm,  2.8  mm,  3.5  mm).  The 
TEG  is  operated  at  the  conditions  of  Cg  =  1,  popt  =  100%, 
qs  -  900  W  m-2,  and  T„  =  25  °C.  The  resistance  of  the  external 
loading  is  constant  in  all  simulations  and  it  is  assumed  to  be  equal 
to  that  of  the  TE  pair  at  32  °C,  that  is,  the  values  of  the  resistance  at 
Geometry  A,  B,  and  C  are  0.0416, 0.0212,  and  0.0116  Q,  respectively. 


3.1.  Comparison  of  real  model  and  equivalent  model 


The  net  output  power  is  obtained  by  subtracting  the  power  con¬ 
sumption  from  the  output  power  of  the  solar  TEG. 

2.6.  Numerical  validation 

By  taking  the  temperature-dependent  material  properties  from 
the  study  of  Li  et  al.  [16],  the  numerical  simulation  was  performed 


Material  properties  and  thicknesses  of  substrate  and  electrode  [27,28], 

Pm-'r1)  y(S  m-1)  Thickness  (mm) 

Substrate  (A1203)  35.3  0  0.635 

Electrode  (Cu)  386  581  x  105  0.06 


The  purpose  of  the  present  study  is  to  develop  an  equivalent 
model  to  simplify  and  accelerate  the  simulation.  Therefore,  the 
predicted  output  power  and  total  efficiency  of  the  TEG  at  Geometry 
B  based  on  the  real  model  are  compared  to  those  based  on  the 
equivalent  model.  In  Fig.  6a,  the  thermal  concentration  ratio  is 
increased  by  enlarging  the  substrate  area.  It  can  be  found  that  the 
equivalent  model  overestimates  the  performance  of  the  TEG.  When 
the  thermal  concentration  ratio  is  low  (i.e.  QjjS40),  the  results  of 
the  equivalent  model  are  close  to  those  of  the  real  model.  The 
difference  between  the  two  models  becomes  obvious  when  the 
thermal  concentration  ratio  is  large  to  a  certain  extent,  even  though 
the  boundary  conditions  between  the  two  models  are  identical.  The 
deviation  between  the  two  models  stems  from  the  thermal 
spreading  resistance  and  the  constriction  thermal  resistance 
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Fig.  5.  Comparisons  of  (a)  the  efficiency  of  TEG  and  (b)  the  output  power  and  effi¬ 
ciency  of  TEG. 


Fig.  6.  Distributions  of  (a)  output  power  and  total  efficiency  of  the  TEG  at  Geometry  B 
and  (b)  relative  error  of  output  power  of  the  thermoelectric  element  sizes. 


caused  in  the  real  model.  When  heat  flows  from  a  small  area  to  a  big 
one,  the  thermal  spreading  resistance  is  induced  [12,30].  Inversely, 
when  heat  moves  from  a  big  area  to  a  small  one,  the  constriction 
thermal  resistance  is  induced  [12,30],  Because  the  substrate  area  in 
the  equivalent  model  is  invariant  when  the  thermal  concentration 
ratio  varies  (Fig.  2a),  the  realistic  behavior  of  thermal  spreading 
resistance  and  constriction  thermal  resistance  are  disregarded.  This 
is  the  reason  why  the  equivalent  model  overestimates  the  output 
power. 

The  profiles  of  the  relative  error  of  output  power  between  the 
two  models  at  the  three  geometries  are  shown  in  Fig.  6b  where  the 
relative  error  is  defined  as 

Realative  error  (%J  fa  Equivalent  -  Eeai  x  10Q  (17) 

The  deviation  is  ranked  by  Geometry  C  >  B  >  A.  This  is  because  the 
area  ratio  between  the  substrate  and  the  electrodes  at  Geometry  C 
is  larger  than  those  at  Geometries  A  and  B.  This  implies,  in  turn,  that 
the  thermal  spreading  resistance  and  constriction  thermal 


resistance  are  more  significant  at  Geometry  C  [30,31],  An  expo¬ 
nential  function,  G  =  exp(mln  Oh  +  n),  is  constructed  to  fit  the 
relative  errors  at  various  geometries  and  the  fitted  curves  are  also 
shown  in  Fig.  6b.  The  values  of  m  and  n  at  the  three  geometries  are 
given  in  Table  2.  Fig.  6b  suggests  that  the  established  exponential 
function  can  precisely  describe  the  relative  error.  Based  on  the 
established  function,  the  predicted  output  power  from  the  equiv¬ 
alent  model  is  modified  by 


(18) 


Following  the  modification,  the  effects  of  the  thermal  spreading 
resistance  and  the  constriction  thermal  resistance  have  intrinsically 


Table  2 

Curve  fitting  functions  (G  =  emln  c,h+")  of  relative  error  between  the  real  model  and 
the  equivalent  model. 

Geometry  ABC 

m  0.656  0.744  1.156 

n  -1.577  -1.173  -1.668 
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Table  3 

A  list  of  substrate  area  and  thermal  concentration  ratio  of  three  geometry  types. 
Substrate  Area  (mm2)  Geometry 


40  x  40 
50  x  50 
60  x  60 
70  x  70 
80  x  80 
90  x  90 
100  x  100 

no  x  no 


50 

78.13 

112.5 

153.13 
200 

253.13 

312.5 

378.13 


25.51 

39.86 

57.40 

78.13 

102.04 

129.15 

159.44 

192.92 


6.38 

9.96 

14.35 

19.53 

25.51 

32.29 

39.86 

48.23 


been  taken  into  account  in  the  equivalent  model.  In  other  words, 
the  disadvantage  of  disregarding  the  thermal  spreading  resistance 
and  constriction  thermal  resistance  in  the  equivalent  model  is 
resolved  from  the  modified  output  power.  After  testing,  the 
established  curves  of  relative  error  can  be  applied  for  the  substrate 
conductivity  ranged  from  35.3  to  130  W  m-1  K-1.  The  following 
simulations  are  performed  based  on  the  equivalent  model  with 
modification. 

3.2.  Effect  of  substrate  area  or  thermal  concentration  ratio 

The  effect  of  substrate  area  or  thermal  concentration  ratio  on 
the  performance  of  the  TEG  is  investigated  in  this  section  where 
eight  different  substrate  areas  (i.e.  the  substrate  types  1—8)  are 
considered,  as  shown  in  Table  3.  The  corresponding  thermal  con¬ 
centration  ratios  of  the  eight  substrate  areas  at  the  three  geome¬ 
tries  are  given  in  the  table  as  well.  The  substrate  area  of 
40  x  40  mm2  is  a  common  size  for  commercial  thermoelectric 


Fig.  8.  Distributions  of  (a)  output  power  and  (b)  total  efficiency  of  the  TEG  at  various 
thermal  concentration  ratios  and  geometries.  (hL  =  1000  W  m~2  K-1). 


Fig.  7.  Distributions  of  (a)  output  power  and  (b)  total  efficiency  of  I 
B)  at  various  thermal  concentration  ratios  and  cooling  conditions. 
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Table  4 

The  values  of  Cth  and  the  corresponding  substrate  areas  of  three  geometry  types. 


Geometry 

A 

B 

C 

Ct„ 

25.51 

39.86 

25.51 

39.86 

25.51 

39.86 

Cross-sectional  area 
of  a  TE  pair  (mm2) 

Substrate  area  (mm2) 

2 

28.57  x  28.57 

35.71  x  35.71 

3.92 

40  x  40 

50  x  50 

15.68 

80  x  80 

100  x  100 

ABC 


Geometry 


Geometry 


Fig.  9.  Profiles  of  (a)  the  output  power  and  output  power  density  and  (b)  the  total 
efficiency  of  the  TEG.  (hL  =  1000  W  nr2  K~’). 


modules  and  it  is  regarded  as  the  smallest  size  among  the  eight 
sizes,  whereas  the  area  of  110  x  110  mm2  (slightly  larger  than  the 
size  demonstrated  in  Fig.  lb)  is  considered  as  the  largest  size.  The 
distributions  of  the  output  power  and  total  efficiency  of  the  solar 
TEG  at  Geometry  B  versus  thermal  concentration  ratio  (or  substrate 
type)  are  shown  in  Fig.  7  where  four  different  cooling  conditions  or 
convection  heat  transfer  coefficients  are  included  for  comparison. 
The  output  power  and  total  efficiency  increase  when  the  thermal 
concentration  ratio  increases  in  all  the  cooling  conditions.  The 
difference  among  the  three  curves  of  hL  =  500,  1000,  and 
2000  W  m  2  K-1  is  almost  imperceptible,  and  they  are  very  close  to 
that  of  the  ideal  case  where  the  convection  heat  transfer  coefficient 
is  assumed  to  be  infinite.  It  follows  that  the  output  power  of  the 
thermal-concentrated  solar  TEG  is  independent  of  the  convection 
heat  transfer  coefficient  when  the  coefficient  is  greater  than 
500  W  itt2  K-1.  For  the  substrate  type  3  (Cth  =  57.40)  at 
hL  =  500  W  m-2  K-1,  if  a  1-D  thermal  resistance  model  is  used  to 
calculate  the  thermal  resistance  of  the  solar  TEG,  the  ratio  of  the 
convectional  thermal  resistance  to  the  total  thermal  resistance  is 
around  2.0%.  The  ratio  decreases  when  hL  increases.  This  reveals 
that  the  convectional  thermal  resistance  is  by  far  smaller  than  the 
total  thermal  resistance.  This  is  the  reason  why  the  output  power  is 
independent  of  the  convection  heat  transfer  coefficient  when  the 
coefficient  is  greater  than  500  W  m-2  I<-1.  The  preceding  charac¬ 
teristic  is  quite  different  from  the  behavior  of  optical-concentrated 
solar  TEG  [16]  which  is  sensitive  to  the  convection  heat  transfer 
coefficient.  This  is  due  to  the  much  larger  substrate  area  of  each  TE 
pair  in  the  investigated  TEG  when  compared  to  that  in  a  commer¬ 
cial  thermoelectric  module.  Therefore,  a  larger  heat  convection  area 
at  the  cold  side  is  obtained  in  the  thermal-concentrated  solar  TEG 
and  heat  can  be  dissipated  from  the  cold  side  to  the  environment 
more  effectively.  When  the  convection  heat  transfer  coefficient  is  as 
low  as  10  W  rrT2  K-1  which  is  at  the  condition  of  free  convection 
[32],  an  obvious  decay  in  the  performance  of  the  TEG  is  observed. 

The  performances  of  the  TEG  at  the  three  geometries  are  pre¬ 
sented  in  Fig.  8  where  the  convection  heat  transfer  coefficient  is 
1000  W  m-2  K-1.  An  increase  in  thermal  concentration  ratio 
monotonically  increases  the  output  power,  regardless  of  the  ge¬ 
ometry  type  (Fig.  8a).  Increasing  thermal  concentration  ratio  also 
monotonically  increases  the  total  efficiencies  of  the  TEG  at  Geom¬ 
etries  B  and  C,  whereas  there  exists  an  optimal  thermal  concen¬ 
tration  ratio  or  substrate  area  (=90  x  90  mm2)  at  Geometry  A 
where  the  maximum  total  efficiency  of  the  TEG  is  4.15%.  This  can  be 
explained  by  the  lower  ZT  value  of  the  thermoelectric  material  at 


A  list  of  properties  of  fins  and  coolants  [32,34], 

Fin  Coolant 

Material  Aluminum  Air 

Property  fcf  =  237  W  m  1  K  1  p1  =  1.1707  kg  m  3 

Ma  =  183.6  x  10-7  N  s  m2 
Cp,a  =  1007Jkg1K-1 


Water 

Pw  =  997.1  kg  m-3 

=  891  x  10~6  N  s  m-2 
Cp,„  =  41 79  J  kg-1  K-1 
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high  temperatures  (Fig.  3d).  The  average  temperature  of  the  ther¬ 
moelectric  element  at  the  substrate  area  of  90  x  90  mm2  from  the 
numerical  prediction  is  about  115  °C  which  is  close  to  the  optimal 
temperature  shown  in  Fig.  3d. 

As  observed  in  Figs.  7  and  8,  the  thermal  concentration  ratio 
(Qh)  is  an  important  parameter  for  the  operation  of  TEG.  Two 
values  of  Cth,  namely,  25.51  and  39.86,  are  considered  for  compar¬ 
ison  and  the  corresponding  substrate  areas  of  the  three  geometries 
are  listed  in  Table  4.  As  shown  in  Fig.  9,  the  TEG  at  Geometry  C  gives 
the  best  performance  among  the  three  geometries  at  Cth  =  25.51 
and  39.86.  At  a  given  Cth,  the  TEG  at  Geometry  C  has  the  larger 
substrate  area  than  those  at  Geometries  A  and  B,  implying  that  a 
larger  input  solar  energy  is  given  in  the  former.  It  is  noteworthy  that 
the  output  power  of  the  TEG  at  Geometry  B  (a  larger  element  cross- 
sectional  area)  is  slightly  larger  than  that  at  Geometry  A,  but  the 
power  density  of  the  former  is  lower  than  that  of  the  latter.  This 
behavior  is  consistent  with  the  results  of  Jang  et  al.  [33],  The  TEGs  at 


Geometries  A  and  B  have  the  same  element  length  (I).  The  results 
reflect  that,  for  a  TEG  at  a  given  element  length,  decreasing  the 
cross-sectional  area  of  the  thermoelectric  element  is  a  feasible 
route  to  improve  the  performance  of  the  TEG,  whether  the  sub¬ 
strate  area  or  Cth  is  given. 

3.3.  Effect  of  cooling  method 

The  theoretical  method  described  earlier  is  used  to  evaluate  the 
net  output  power  of  the  solar  TEG,  and  two  different  coolants,  air 
and  water,  are  considered.  The  properties  of  the  fins  and  coolants 
are  listed  in  Table  5  [35],  It  can  be  seen  that  fin-to-fin  spacing  (Dg) 
has  the  highest  order  in  Eq.  (16),  so  it  is  chosen  as  the  operating 
parameter.  The  height  (Hf)  and  thickness  (tf)  of  the  fins  are  14  mm 
and  0.2  mm,  respectively  [24],  The  base  area  of  the  heat  sink  is 
assumed  to  be  equal  to  the  substrate  area  of  the  thermoelectric 
module,  so  the  width  (WHS)  and  the  length  (IHs)  of  the  sink  are 
decided  by  the  given  substrate  area.  The  distribution  of  the  power 
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Geometry  B 
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-12o 


1.6  2  2.4  2.8  3.2 


Dg  (mm) 


Dg  (mm) 


Fig.  10.  Distributions  of  (a)  pumping  power  for  water  cooling  and  (b)  the  power  of  the  Fig.  11.  Distributions  of  net  ou 

TEG  at  Geometry  B.  57.40  and  (b)  102.04. 
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for  the  TEG  at  Geometry  B  is  shown  in  Fig.  10  where  the  water 
cooling  is  carried  out.  Fig.  10a  shows  that  the  pumping  power  de¬ 
creases  when  Dg  increases  or  Cth  decreases.  Fig.  10b  suggests  that 
the  required  pumping  power  is  substantially  smaller  than  that  of 
the  output  power  of  the  TEG  within  the  investigated  range  of  Dg 
(Cth  —  159.44).  Consequently,  the  net  output  power  is  almost 
invariant  with  varying  Dg  when  water  cooling  is  adopted.  The 
comparisons  of  net  output  power  between  water  cooling  and  air 
cooling  at  Cth  =  57.40  and  102.04  are  shown  in  Fig.  11.  With  the 
condition  of  Qh  =  57.40  and  under  the  same  cooling  condition 
(hi  —  1000  W  m-2  K_1),  a  negative  net  output  power  is  obtained  for 
air  cooling  when  a  small  Dg  is  used.  The  positive  net  output  power 
for  air  cooling  can  be  obtained  only  when  Dg  is  large  enough 
(Fig.  11a).  However,  when  Cth  is  102.04,  the  positive  net  output 
power  under  air  cooling  cannot  be  obtained  in  the  investigated 
range  of  Dg  (Fig.  lib).  Additionally,  the  Reynolds  number  will 
exceed  the  laminar  limit  (2300)  [35]  if  Dg  is  too  large  (the  limit 
value  of  Dg  becomes  smaller  when  Cth  increases  as  shown  in 
Fig.  lib),  which  doesn’t  take  place  for  water  cooling.  Accordingly, 
water  cooling  is  a  better  cooling  method  to  maximize  the  net 
output  power  of  the  solar  TEG,  stemming  from  the  larger  specific 
heat  of  water. 

4.  Conclusions 

The  performance  of  a  thermal-concentrated  solar  TEG  has  been 
investigated  using  a  numerical  method  where  the  temperature- 
dependent  properties  of  the  commercial  thermoelectric  material 
were  taken  into  account.  In  the  method,  an  equivalent  model  was 
developed  to  simulate  the  real-sized  TEG  module  and  an  analytical 
method  was  used  to  evaluate  the  power  consumption  for  cooling 
the  TEG  module.  By  introducing  an  exponential  function  with  the 
consideration  of  the  constriction  thermal  resistance  and  thermal 
spreading  resistance  in  the  equivalent  model,  numerical  simula¬ 
tions  with  much  less  computational  time  was  achieved. 

Three  different  geometry  types  (Geometries  A,  B,  and  C)  of  TEG 
were  tested  and  compared  with  each  other.  The  predictions  suggest 
that  the  output  power  of  the  solar  TEG  increases  with  increasing 
substrate  area,  stemming  from  the  increased  thermal  concentration 
ratio.  The  smallest  element  (Geometry  A)  gives  the  best  perfor¬ 
mance  and  the  maximum  efficiency  is  4.15%  when  the  substrate 
area  is  90  x  90  mm2,  as  a  consequence  of  lower  ZT  value  of  the 
thermoelectric  materials  at  high  temperatures.  When  the  element 
length  is  fixed,  decreasing  the  cross-sectional  area  of  the  thermo¬ 
electric  element  is  a  better  route  to  improve  the  performance.  The 
cooling  conditions  almost  play  no  part  in  the  performance  of  the 
solar  TEG  when  the  forced  convection  with  the  conditions  of 
hL  a  500  W  m-2  K-1  is  encountered.  According  to  the  theoretical 
method,  a  larger  and  positive  net  output  power  can  be  achieved 
using  water  cooling  for  the  solar  TEG.  This  work  has  provided  a 
useful  insight  into  the  design  of  practical  thermal-concentrated 
solar  TEG  system. 
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